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ABSTRACT 



A computer program (written in FORTRAN) is presented 
which uses Monte Carlo techniques to simulate one-searcher, 
one-target passive acoustic ASW search that terminates at 
detection. A threshold crossing detection model is used, 
and stochastic variations in the acoustic signal are modeled 
using either a Lambda-Sigma Jump or Gauss-Markov error 
process. Both platforms have the capability of detecting 
each other, and area and barrier searches are modeled. 
Features of the program include interactive data input, 
extensive use of graphical displays , and thorough statis- 
tical analysis of the results of the simulation. 
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I. 



DESCRIPTION OF THE SIMULATION MODEL 



A. INTRODUCTION 

Although the body of analysis devoted to search theory 
is extensive (see, for instance, the lists of references in 
Koopman [Ref. 1] and Washburn [Ref. 2]), most of the tech- 
niques are applicable only to stationary targets. Recently 
there have been some important advances in the theory of 
search for moving targets (see Eagle [Ref. 3] , Brown [Ref. 4] , 
Stewart [Ref. 5] , Stone [Ref. 6] , and Washburn [Ref. 7]) . 

Even so, it is still a very difficult, if not impossible, 
computational problem to calculate the probability of 
detection or the mean time to detection when the searcher 
has a speed close to that of the target and follows a 
realistic (i.e. not random) track. As a result, computer 
simulation of the search process is often used to evaluate 
the effectiveness of search tactics. 

The Monte Carlo simulation computer program presented in 
this thesis is an attempt to model the passive acoustic 
search process in a one-searcher, one-target scenario where 
both the searcher and target platforms have sensors capable 
of detecting each other. The results of the simulation 
allows the analyst to closely approximate parameters (e.g. 
cumulative probability of detection, detection times, detec- 
tion ranges) which are essential in planning search 
operations. With these results, the military planner can 
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make informed decisions concerning asset allocation and 
specific platform employment. 

The computer program, Passive Acoustic Search Simulation 
(PASS) , is based on an acoustic fluctuation model described 
by McCabe [Ref. 10:pp.7-23] and Hurley [Ref. 11 : pp . 25-34 ] . 
Precedessor programs include the FORTRAN program BEAR, and 
the ALGOL program SEARCH. BEAR was written by L. K. Arnold 
of D. H. Wagner, Associates (see Arnold [Ref. 8] for program 
documentation). SEARCH, a derivitive of BEAR, was written 
by W. J. Browning and J. Risberg, also of Wagner, 

Associates (see Browning [Ref. 9] for program documentation). 
SEARCH has been used extensively in operational and exercise 
environments under the direction of Submarine Development 
Squadron Twelve (SUBDEVRON-12 ) , where it is currently 
resident on the Semi-Automatic Reconstruction Facility 
(SARF) computer. The basic structure of PASS is similar to 
that of BEAR, with some of the features of SEARCH added. 
Additionally, some modifications were made to algorithms in 
these earlier programs, and new features were added, 
including interactive data input and extensive use of 
graphical presentation of simulation results . 

PASS is written in FORTRAN and is designed for use at 
the Naval Postgraduate School (NPGS) . PASS is not portable 
in that it uses external subroutines in the Non-International 
Mathematics and Statistics Library (NONIMSL) and the DISSPLA 
graphics system. The subsequent sections of this chapter 
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are designed to give the reader an overview of the model 
without going into the specifics of the program code or 
logic. If the reader is interested only in learning how to 
use PASS, it is recommended that this chapter and Appendix A 
be studied. Appendix A is designed as a User's Guide and 
does not require any knowledge of the program code other 
than variable names, which are described in Appendix B. If 
the reader is interested in the structure of the program and 
the model, then this chapter and Chapter II should be studied. 
Some examples of varied applications of PASS are shown in 
Chapter III. 

B. GENERAL PROGRAM DESCRIPTION 

Program PASS is a Monte Carlo simulation of two passive 
sonar platforms in a search scenario. Arbitrarily, one 
platform is called the "searcher", and one the "target", 
where the target is typically thought of as a submarine, and 
the searcher as any passive sonar platform. Since both ships 
have passive sensors, the assignment of the names searcher 
and target does not connotate strict roles (i.e., the searcher 
may be detected first, and becomes a target in the truest 
sense I ) . 

During each replication, the target's motion is 
partially random within a user-defined search area. The 
searcher motion is deterministic and controlled by user 
input. Each sensor is subject to random acoustic fluctuations 
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which induce random variations in detection ranges over 
time. Detection occurs whenever one platform obtains a 
positive signal excess, unless the user employs an integra- 
tion time model, which requires weak signals to be present 
longer than strong signals for detection to occur. 

The output of the program consists of distributions of 
parameters of interest (e.g. time of detection and range of 
detection) , and point estimates derived from the statistics 
of these distributions (e.g. mean time to detection and 
mean detection range) . The format of the output is numerical 
for point estimates, and graphical for distributions. 

C. SEARCH SCENARIOS USED IN PASS 

1 . Area Search 

The target is constrained to move within a specific 
rectangular area ("search area") in a "semi-random" fashion. 
Target speeds are distributed uniformly over a user speci- 
fied speed range, and courses are uniform on (0, 360) 
degrees. Times between course and speed changes are exponen- 
tially distributed. The size of the search area and the 
parameters of the random motion are user specified. The 
searcher moves along a track which need not be contained 
within the search area. The searcher track and speeds are 
user specified. 

2 . Barrier Search 

The target is constrained laterally to be within 
certain user specified bounds, thus simulating a "choke point". 
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Target motion through the choke point is generally vertical 
(i.e. from "top" to "bottom") with random course and speed 
deviations parameterized by the user. The searcher track 
is specified as in Area Search. 

The distribution of the initial lateral position of 
the target can be specified by the user as either uniform 
across the choke point, or with specific probabilities as a 
function of distance across the choke point. In using the 
latter option, and by limiting the maximum course deviation, 
the user can control the distribution of the lateral target 
position when the target penetrates the barrier. 

D. THE ENVIRONMENTAL MODEL 

The environment is modeled in a conventional manner. 

PASS requires range/propagation loss information for both 
the searcher and target sensors. Direct path (DP) and 
convergence zone (CZ) data are entered separately. A cubic 
interpolation routine is used to determine propagation loss 
between the user-input data points. The convergence zones 
are modeled as inverted "square-wells" superimposed over the 
DP propagation loss curve. Up to five convergence zones are 
allowed. 

E. THE DETECTION MODEL 

In PASS a detection is always a "secure detection" . That 
is, a detection by the searcher means that the searcher 
detected the target before the target detected the searcher. 
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Simultaneous detections are treated as special cases where 
neither the searcher or the target gets credit for a 
detection. The simulation terminates (replication ends) 
whenever: 

1. The searcher makes a detection. 

2. The target makes a detection (sometimes referred to 
as a counter-detection) . 

3. A simultaneous detection occurs. 

4. In the barrier scenario, the target crosses a user 
specified lower boundary. 

5. User specified maximum search time is exceeded. 

The propagation loss curves are used in conjunction with 
searcher and target f igure-of-merit (FOM) to determine 
detection ranges . The program uses FOM data which must be 
precalculated by the user based on specific platform source 
levels (Ls) , backgroung noise (self noise and ambient noise) 
and sensor directivity (Le) , and processor (machine and 
operator) recognition differential (Nrd) . 

1 . Searcher Figure-of-Merit 

The Ls of the target is assumed to be independent 
of target speed. In effect flow related noise and speed 
dependent propulsion noise are ignored. At the user's 
option, target Ls may be made a function of target aspect. 

It is assumed that the deterministic sensor performance of 
the searcher is affected by the searcher speed. This results 
in a searcher FOM that is dependent on searcher speed, and, 
optionally, target aspect. 
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2 . Target Figure-of-Merit 



The Ls of the searcher is assumed to be speed depen- 
dent, and independent of searcher aspect. This results in a 
target FOM that is dependent on searcher speed. 

3 . The Threshold Crossing Model With a Stochastic Error 

Process 

Given a FOM against a target for a search platform, 
the maximum detection range of the searcher is determined 
from the propagation loss curves . This results in the sensor 
conforming to the "definite range law", or "cookie-cutter" 
performance. This type of sensor readily lends itself to 
mathematical analysis, but has the disadvantage of not 
reflecting actual sensor performance. What is seen under 
operational conditions is a variation in signal strength 
(and thus a variation in detection ranges) over time, even 
for targets at fixed range and source level. 

To model these variations in received acoustic 
signal, a stochastic error process is added to the determin- 
istic figure-of-merit . This results in the signal excess 
consisting of deterministic and random components. For a 
detailed discussion of threshold crossing models and 
associated stochastic error processes, see McCabe [Ref. 10] 
and Hurley [Ref. 11] . See Appendix D for a discussion of 
the relationship of signal excess (SE) to f igure-of-raerit 
(FOM) . 

The mean signal excess (SE) is defined as the 
expected difference between the signal-to-noise ratio (SNR) 
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at the hydrophone output (in decibels, db) and the SNR that 
is estimated to afford a probability of detection of 0.5. 

The model used for the "observed" signal excess, as 
a function of time in PASS, is of the form: 

SE = SE + X + Y (1.1) 

where SE is the deterministic mean signal excess, and X and 
Y are random variables that describe the fluctuations in the 
signal level. 

If we let the subscript s apply to the searcher and 
t for target, the model used in PASS is: 

SE g (t) = SEg ( t) + X(t) + Y (t) (1.2) 

SE t (t) = SE t (t) + Y ( t) + Z(t) (1.3) 

where : 

1. X(t) is a stochastic error process which describes the 
fluctuations "local" to the searcher. 

2. Z(t) is a stochastic error process which describes the 
fluctuations "local" to the target. 

3. Y(t) is a stochastic error process which describes 
acoustic fluctuations of a global nature whose effect 
is common to the searcher and target. 

4. X ( t ) , Y ( t ) , and Z(t) are independent processes. 

X(t) and Z(t) may be thought of as modeling the 
fluctuations that affect only the performance of the 
detecting platform, such as onboard fluctuations in sensor/ 
processor performance, including operator performance. 
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Y(t) may be thought of as modeling the gross 
acoustic fluctuations that occur along the "environmental 
path" between the platforms, such as temperature and 
salinity gradients and changes in bottom conditions. 

With signal excess now a function of the mean signal 
excess and the stochastic effor processes, we have the means 
to model the behavior of the sensor by employing a threshold 
crossing detection model as follows: 

1. If SE(t) is less than 0, the probability of detection 
is 0.0. 

2. If SE(t) is greater than or equal to 0, the probability 
of detection is 1.0. 

Thus, a detection occurs at time t if, and only if, SE(t) is 
greater than or equal to zero. The fact that SE(t) is now 
undergoing random fluctuations about the mean signal excess 
results in random variations in the sensor detection range. 

F. STOCHASTIC ERROR PROCESSES 

PASS allows the user to choose either a Lambda-Sigma 
Jump (LSJ) or a Gauss-Markov (GMA) error process to model 
the acoustic fluctuations. Hurley [Ref. 11] gives an over- 
view of these processes, and others used in search simulation 
programs. McCabe [Ref. 10] provides a more mathematically 
complete treatment of the processes, along with a compre- 
hensive comparison of the LSJ and GMA models. In PASS we 
use the sum of two LSJ or GMA process, frequently referred 
to as a compound process. In the following discussion, a 
simple process, vice the compound process is considered for 
clarity . 
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1 . The Lambda-Sigma Jump Error Process 



Define a Lambda-Sigma Jump error process, £(t), as 

follows : 



£ (t+s) 
(six) 

z (t) = 



= z (t) £ ( t) + [l-z (t) ] n 




if s 



< T 



0 if s = x 

V- 



(1.4a) 



(1.4b) 



where : 

1. x is an exponential random variable with rate parameter X. 

2. r) is a normal random variable with zero mean and 

variance cr 2 

The process can be though of as a "marked" Poisson 
process with rate parameter X, where the magnitude of the 
marks are themselves random variables with a normal distri- 
bution (with zero mean and variance a 2 ) . The standard 
deviation of the normal distribution, a, can be thought of 
as a scale parameter in this process. The figure-of-merit , 
then, is constant over exponentially spaced time intervals, 
and changes to new levels determined by a normal distribution 
at the end of each of these time intervals . 

The covariance function for £(t), is: 

Cov [£(t), 5 (t+s) ] = a 2 e" Xs (1.5) 



This function is a measure of the correlation between values 
of £(t) at different times. From this we can see that £(t) 
is second order stationary. That is, the covariance depends 
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on the time difference, s, and not the time, t. The para- 
meter X determines the amount of dependence between successive 
acoustic levels, and can be chosen to yield any result from 
independence to complete dependence. The resulting process 
is then a piecewise continuous function in which there are 
periods of complete dependence of detection opportunity 
interrupted by fluctuations that introduce independence 
between these periods. 

In PASS, then, we must provide a total of six param- 
eters to completely describe the compound error processes 
used to model the acoustic fluctuations. They are: 

1. X x = rate parameter for the searcher-local error 
process , X ( t) . 

2. = scale parameter for the searcher-local error 
process , X ( t) . 

3. X 2 = rate parameter for the global error process, Y(t). 

4. <J 2 = scale parameter for the global error process, Y(t). 

5. X 3 = rate parameter for the target local error 
process , Z ( t) . 

6. a 3 = scale parameter for the target local error 
process , Z ( t) . 

The selection of appropriate values for X and a is a 
not well understood function of the environment, sensor, and 
processor. In practice, the selection of the values are 
subjective decisions based on the experience of the analysist 
as what seems to yield reasonable results. Tehan [Ref. 12] 
provides a discussion of estimation techniques for 
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